Gravitational Lensing by Elliptical Galaxies, and 
the Schwarz Function 

CD. Fassnacht, C. R. Keeton and D. Khavinson 



Abstract. We discuss gravitational lensing by elliptical galaxies with some 
particular mass distributions. Using simple techniques from the theory of 
quadrature domains and the Schwarz function (cf. [18]) we show that when 
the mass density is constant on confocal ellipses, the total number of lensed 
images of a point source cannot exceed 5 (4 bright images and 1 dim image). 
Also, using the Divo-Nikliborc converse of the celebrated Newton's theorem 
concerning the potentials of ellipsoids, we show that "Einstein rings" must 
always be either circles (in the absence of a tidal shear), or ellipses. 



1. Basics of gravitational lensing 

Imagine n co-planar point-masses (e.g., condensed galaxies, stars, black holes) that 
lie in one plane, the lens plane. Consider a point light source S (a star, a quasar, 
etc.) in a plane (a source plane) parallel to the lens plane and perpendicular to the 
line of sight from the observer, so that the lens plane is between the observer and 
the light source. Due to deflection of light by masses multiple images 5i, 5*2, . . . of 
the source may form (cf. Fig. 1). Fig. 2 and Fig. 3 illustrate some further aspects 
of the lensing phenomenon. 
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Figure 2. Lcnsing of a galaxy by a cluster of galaxies; the blue 
spots are all images of a single galaxy located behind the huge 
cluster of galaxies. (Credit: NASA, W. N. Cooley (Princeton), E. 
Turrer (Princeton) and J. A. Tyson (AT&T and Bell Labs).) 
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Figure 3. The bluish bright spots are the lensed images of a 
quasar (i.e., a quasi-stellar object) behind a bright galaxy in the 
center. There are 5 images (4 bright + 1 dim), but one cannot 
really see the dim image in this figure. (Credit: ESA, NASA, K. 
Sharon (Tel Aviv University) and E. Ofek (Caltech).) 



2. Lens equation 

In this section we are still assuming that our lens consists of n point masses. 
Suppose that the light source is located in the position w (a complex number) in 
the source plane. Then, the lensed image is located at z in the source plane while 
the masses of the lens L are located at the positions z^, j = 1, . . . , n in the lens 
plane. The following simple equation, obtained by combining Fermat's Principle 
of Geometric Optics together with basic equations of General Relativity, connects 
then the positions of the lensed images, the source and the positions of the masses 
which cause the lensing effect 



(2.1) 



where cTj are some physical (real) constants. For more details on the derivation 
and history of the lensing equation ( |2.1[ ) we refer the reader to [19], [12], [14], [21]. 
Sometimes, to include the effect caused by an extra ("tidal") gravitational pull by 
an object (such as a galaxy) far away from the lens masses, the right-hand side of 



(2.1) includes an extra linear term 72, thus becoming 

n 

- E - rz, (2.2) 



w 
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where 7 is a real constant. The right-hand side of (2.1 1 or (2.2) is called the lensing 
map. The number of solutions z of (2.1 1 (or (2.2 1) is precisely the number of images 

n 

of the source w generated by the lens L. Letting r{z) = 'Y^ + 72 + iS, the lens 

1 ^ 

equations (2.1) and (2.2) become 

z — t(z) — 0, 

where r(z) is a rational function with poles at Zj, j = 1, . 



(2.3) 

, , n and infinity if 7 7^ 0. 



3. Historical remarks 

The first calculations of the deflection angle by a point mass lens, based on New- 
ton's corpuscular theory of light and the Law of Gravity, go back to H. Cavendish 
and J. Michel (1784), and P. Laplace (1796)— cf. [20]. J. Soldner (1804) — cf. [21] 
is usually credited with the first published calculations of the deflection angle and, 
accordingly, with that of the lensing effect. Since Soldner's calculations were based 
on Newtonian mechanics they were off by a factor of 2. A. Einstein is usually given 
credit for calculating the lensing effect in the case of n = 1 (one mass lens) around 
1933. Yet, some evidence has surfaced recently that he did some of these calcula- 
tions earlier, around 1912 — cf. [17] and references therein. The recent outburst 
of activity in the area of lensing is often attributed to dramatic improvements of 
optics technology that make it possible to check many calculations and predictions 
by direct visualization. 

H. Witt [24] showed by a direct calculation that for n > 1 the maximum 
number of observed images is < -I- 1. Note that this estimate can also be derived 
from the well known Bezout theorem in algebraic geometry (cf. [3,8,9,22]). In [11] 
S. Mao and A. O. Fetters and H. J. Witt showed that the maximum possible 
number of images produced by an n-lens is at least 3n + 1. A. O. Fetters in [13], 
using Morse's theory, obtained a number of estimates for the number of images 
produced by a non-planar lens. S. H. Rhie [15] conjectured that the upper-bound 
for the number of lensed images for an n-lens is 5n — 5. Moreover, she showed 
in [16] that this bound is attained for every n > 1 and, hence, is sharp. Rhie's 
conjecture was proved in full in [8]. Namely, we have the following result. 

Theorem 3.1. The number of lensed images by an n-mass, n > I, planar lens 
cannot exceed bn — b and this bound is sharp [16]. Moreover, the number of images 
is an even number when n is odd and odd where n is even. 

The proof of the above result rests on some simple ideas from complex dy- 
namics (cf. [9, 10]). 



4. "Thin" lenses with continuous mass distributions 

If we to replace point masses by a general, real-valued mass distribution /i, a 
compactly supported Borel measure in the lens plane, the lens equation with shear 
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(2.2) becomes 



/ ^^^(0 - (A ^\ 

w = z— I -—jz. (4.1) 



Here is a bounded domain containing the support of 



11. T he case of the atomic 

n 

measure fi — J^'^j^zjj fJi e M is covered by Theorem 3.1 Also, as noted in [8], 

1 I I 

if we replace n-point-masses by n non-overlapping radially symmetric masses, the 

total number of images outside of the region occupied by n-masses is still 5n — 5 

when 7 = 0, and < 5n when ^ =/= 0. The reason for that, of course, is that the 

Cauchy integral 

'^^(^^ \z-z,\>R 



z-C 

for any radially symmetric measure ^ — ^{\C, ~ Zj\) is immediately calculated to 
be equal ^3^, where c is the total mass of the disk : |C ^ ^il < -R}; hence 
reducing this new situation to the one treated in Theorem |3.1[ 

Here is another situation that can be treated with help from Theorem |3.1| 
Recall that a simply-connected domain is called a quadrature domain (of 
order n) if is obtained from the unit disk D := {z : \z\ < 1} via a conformal 
mapping that is a rational function of degree n, = 1^9(15). Of course, all poles 
(3j, j — I, . . . ,n of if will lie outside D. Then if, say, fi is a, uniform mass distribution 
in f2, i.e., = const dxdy, the Cauchy potential term in (4.1 1 for z ^ becomes 

where the coefficients Cj are determined by the quadrature formula associated with 
n (cf. [18] for details). 

Hence, substituting (4.2 1 into (4.11 we again obtain that for such thin lens 
f2 with a uniform density distribution, the number of "bright" images outside f2 
cannot exceed 5n — 5 when no shear is present, or 5n otherwise. 

In this general context the only previously known (to the best of our knowl- 
edge) result is the celebrated Burke's theorem [2] 

Theorem 4.1. A (finite) number of images produced by a smooth mass distribution 
fi is always odd, provided that 7 = (no shear). 

An elegant complex-analytic proof of Burke's theorem can be found in [19]. 
The crux of the argument is this. Take w = and let n+, n_ denote, respectively, 
the number of sense-preserving and sense-reversing zeros of the lens map in (4.1) 
(7-O). 

The argument principle applies to harmonic complex- valued functions in the 
same way it does to analytic functions. Since the right-hand side of (4.1 1 behaves 
like 0{z) near 00, the argument principle then yields that 1 = — n_. Thus, 
giving us the total number of zeros N := + n_ = 2n_ + 1, an odd number. 



6 



Fassnacht, Keeton and Khavinson 



5. Ellipsoidal lens 

Suppose the lens :— |f2+f2<l,a>^>o| is an ellipse. First assume the 
mass density to be constant, say 1. Let c : — — he the focal distance of O. 



The lens equation (4.11 can be rewritten as 

^__WcM(C)_ 

TT J Z-Q 



7z = w, 



(5.1) 



where dA denotes the area measure. Using complex Green's formula (cf. e.g., [19]), 



we can rewrite (5.1 1 for z € C \ as follows: 

1 f (dC 

c 



27ri J z 

on 



7z 



(5.2) 



As is well-known [18], the (analytic) Schwarz function S{C) for the ellipse defined 
by 5(0 = ( on dil can be easily calculated and equals 

-2 , 2ab 



2ab 



(5.3) 



5l(C) + ^2(C)- 



Note that is analytic in il, while S2 is analytic outside H. and 6*2(00) = 0. This 
is, of course, nothing else but the Plemelj-Sokhotsky decomposition of the Schwarz 



function S{Q of dfl. From (5.3) and Cauchy's theorem we easily deduce that for 
z e C \ the lens equation (5.2 1 reduces to 



(5.4) 



Squaring and simplifying, we arrive from (5.4) at a complex quadratic equation 

2ab \ 2a2&2 



z — 7 I zw 



which is equivalent to a system of two irreducible real quadratic equations. Bezout's 
theorem (cf. [8,9], [10], [3]) then implies that (5.1 1 may only have 4 solutions z ^ il. 



For z G n, using Green's formula and (5.3 1 we can rewrite the area integral in (5.1 ) 

dAiC) 



^z + 



1 

27ri 



-z + 



1 



2TTi J C 

on 



[51(0+^2(0] 



-z + S'i(z) = — z + 



2^ + b^ 



(5.5) 



2ab 
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We have used here that the Cauchy transform of S2 |an vanishes in H. since S2 is 
analytic in C \ and vanishes at infinity. Substituting (5.51 into (5.11, we arrive 
at a linear equation 



- 2ab 



(5.6) 



for z £ n. Equation (5.6 1, of course, may only have one root in ft. Thus, we have 
proved the following 



Theorem 5.1. An elliptic lens f2 {say, a galaxy) with a uniform mass density may 
produce at most four "bright" lensing images of a point light source outside and 
one ("dim") image inside f2, i.e., at most 5 lensing images altogether. 



This type of result has actually been observed experimentally - cf. Fig. 4, 
where four bright images are clearly present. It is conceivable that the dim image 
is also there but we can't see it because it is perhaps too faint compared with the 
galaxy. Of course, one has to accept Fig. 4 with a grain of salt since we do not 
expect "real" galaxies to have uniform densities. A model of an elliptical lens, with 
shear, that produces five images (4 bright + Idim) is given in Fig. 5. 





r 






Gravitational Lenses 

pac94-43 ■ ST Sa 0«3 . Oeisw IS. K. 


HST ■ WFPC2 

R.1|n^lUnD»lJHU|.r?AS* 



Figure 4. Four images of a light source behind the elliptical 
galaxy. (Credit: NASA, Kavan Ratnatunga (Johns Hopkins Uni- 
versity) . ) 
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We can extend the previous theorem for a larger class of mass densities. 



Denote by q{x, y) ^ 



1 the equation of T dil. Let q\{x, y) — 



— 1, —b^ < A < stand for the equation of the boundary Vx := of the 
ellipse r^A confocal with fi. 

The celebrated MacLaurin's theorem (cf. [7]) yields that for any z e C \ f2 



1 



dA{0 



1 



Ajcea{Q.x) J ( — z Area(ri) J C~z 



dAiO 



(5.7) 



Thus, if we denote by u{z, A) the Cauchy potential of il\ evaluated at z e C | 

u{z,X) = c{X)un{z,0), 

Areaiflx) _ (a' + A)'/' {b^ 



we obtain from (5.7 1 
where 



c(A) = 



A) 



1/2 



Hence, 



Area(ri) 

duxiz,X) 
dX 



ab 



c'(A)uo(z). 



(5.8) 
(5.9) 



(5.10) 



So, if the mass density /i(A) in only depends on the elliptic coordinate A, i. e., 
is constant on ellipses confocal with inside il, its potential outside equals 



Uf,,n{z) = cun{z). 



The constant c is easily calculated from (5.9)-(5.10l and equals 





c== / fi{X)c\X)dX. 



(5.11) 



(5.12) 



-62 
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It is, of course, natural for physical reasons to assume that /i(A) "f cx3 at the "core" 
of (i.e., when A i —b^), the focal segment [— c, c]. Yet, from (5.121 since (5.9) 

yields c'(A) = O ^(6^ + A) ^^^^ near Aq = —6^, it follows that ^(A) should not 

^1/2 



diverge at the core faster than say 0(^(6^ + A) "'''"j for some positive e, so 

the integral (5.121 converges. Substituting ( 5.11[ ) into the lens equation (5.1 1 with 
constant density replaced by the density fJ.{X) and following again the steps in 



(5.2)-(5.4l we arrive at the following corollary. 



Corollary 5.1. An elliptic lens il with mass density that is constant inside on 
the ellipses confocal with may produce at most four "bright" lensing images of 
a point light source outside fl. 



6. Einstein rings 



For a one-point mass at zi lens with the source at w = the lens equation (2.1 ) 
without shear becomes 

z-^-^ = 0. (6.1) 

z — Zl 



As was already noted by Einstein (cf. [12,19,21] and references cited therein), (6.1 ) 
may have two solutions (images) when zi ^ and a whole circle ("Einstein ring") 
of solutions when zi = 0, in other words when the light source, the lens and the 
observer coalesce - cf. Fig. 6 and Fig. 7. 



Figure 6. A model of an elliptical Einstein ring surrounding an 
elliptical lens with axis ratio 0.5 and uniform density 2. The shear 
in this case must be specially chosen to produce the ring instead 
of point images. Note that the ring is an ellipse confocal with the 
lens - cf. Thm. 6.1 . 
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J095629. 








Figure 7. Einstein rings. The sources in these observed "realis- 
tic" lenses are actually extended, and that is why we see some- 
times arcs rather than whole rings. (Credit: ESA, NASA and 
the SLACS survey team: A. Bolton (Harvard / Smithsonian), S. 
Buries (MIT), L. Koopmans (Kapteyn), T. Treu (UCSB), and L. 
Moustakas (JPL/Caltech).) 



As the following simple theorem shows the "ideal" Einstein rings are limited 
to ellipses and circles in much more general circumstances. 

Theorem 6.1. Let f2 be any planar ("thin") lens with mass distribution /ie. // 
lensing of a point source produces a bounded "image" curve outside of the lens Q., 
it must either be a circle when the external shear 7 = or an ellipse. 

Proof. First consider a simpler case when 7 = 0. If the lens produces an image 



curve r outside il, the lens equation (4.1 1 becomes 



z-w = j (6-2) 

a 

for all z € r. Note that F being bounded and also being a level curve of a harmonic 
function must contain a closed loop surrounding Q. [22]. Without loss of generality. 



we still denote that loop by F. The right-hand side f{z) of (6.2 1 is a bounded 
analytic function in the unbounded complement component of F that vanishes 
at infinity. Hence {z — w)f{z) is still a bounded and analytic function in equal 
to \z — > on F := dQoo- Hence {z — w)f{z) = const and F must be a circle 
centered at w. 
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Now suppose 7 7^ 0. Once again we shall still denote by F a closed Jordan 
loop surrounding il. Denote by the interior of F, f2oo = C \ clos (n^ . Also, by 
translating we can assume that the position of the source w is at the origin. 



The equation (4.1 1 now reads 



dKC) 



jz, z eT. 



(6.3) 



In other words the right-hand side of (6.3) represents the Schwarz function 5 of F, 
analytic in C \ supp/i with a simple pole at oo. It is well-known that this already 
implies that F must be an ellipse (cf. [4,18]) and references therein. For the reader's 
convenience we supply a simple proof. 



Applying Green's formula to ( |6.3| yields (cf. (5.1 l-(5.2l) that for all z G il, 
rfM(C) ^ 1 fdAiO 

Z-( TT J Z-( 



(6.4) 



Let 



h{z) := 



1 



dA{C) 



z eQ. 



(6.5) 



Then, h(z) is analytic in (cf. (5.2 1 ) and, in view of (6.3) and (6.4) 



2|p = 72; |r 



(6.6) 



Thus, h(z) is a linear function and since (6.5) for z e f2 

1 



Hz) 



■ grad 



^ J log|z-C|dA(C)-|z|^ 



we conclude from (6.7) that the potential of 

1 

2^ 



Ufiiz) 



log\z-C\dAiO, zen 



(6.7) 



equals to a quadratic polynomial inside il. The converse of the celebrated theorem 
of Newton due to P. Dive and N. Nikliborc (cf. [7, Ch. 13-14] and references 
therein) now yields that must be an interior of an ellipse, hence F dCl = dfloo 
is an ellipse. □ 

Remark 6.1. One immediately observes that since the converse to Newton's theo- 
rem holds in all dimensions the last theorem at once extends to higher dimensions 
if one replaces the words "image curve" by "image surface" . 
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7. Final remarks 



1 . The densities considered in §5 are less important from the physical viewpoint 
than so-called "isothermal density" which is obtained by projecting onto the 
lens plane the "realistic" three-dimensional density ~ 1/p^, where p is the 
(three-dimensional ) distance from the origin. This two-dimensional density 
could be included into the whole class of densities that are constant on all 
ellipses homothetic rather than confocal with the given one. The reason for the 
term "isothermal" is that when a three-dimensional galaxy has density ~ 1 / p'^ 
the gas in the galaxy has constant temperature (cf. [5] and the references 
therein). 

Recall that the Cauchy potential of the ellipse il:— <.^ + j^<l,a>b>0 



outside of (cf. (|5.2|-(5.4|) equals 



(7.1) 



where k — 2ah/c^ is a constant. Replacing the ellipse by a homothetic 
ellipse nt-.^tn-.^ ^ + ^ <t^, < t < 1. We obtain using ([tT]) for z ^Vl: 

dA{0 . fdAiC) 



u{z, t) 



tn 



— tu 



^ ' t(- z 



Thus, 



du{z, t) 
dt 



,2^2 



(7.2) 



(7.3) 



So, if the "isothermal" density n = j an dilt inside (ignoring con 



stants), we get from ( 7.3 1 that the Cauchy potential of such mass distribution 
outside equals 



u^{z) := Co 



dt 

Vz^ - cH^ ' 



(7.4) 



where the constant Co depends on only. This is a transcendental function 
(one of the branches of arcsin-), dramatically different from the algebraic 



potential in (7.1 1. The lens equation (4.1 )now becomes 





(7.5) 



To the best of our knowledge the precise bound on the maximal possible 



number of solutions (images) of ( 7.5 ) is not known. Up to today, no more than 
5 images (4 bright -1-1 dim) have been observed. However, in [5] there have 
been constructed explicit models depending on parameters a, b and < 7 < 1 



having 9 (i. e., 8-1-1) images. The equation (7.5 1 essentially differs from all 
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the lens equations considered in this paper since it involves estimating the 
number of zeros of a transcendental harmonic function with a simple pole at 
oo. At this point, we are even reluctant to make a conjecture regarding what 
this maximal number might be. 

Note, that in case of a circle ft = {x^ + < l} with any radial density 
H := (p{r), r = x"^ + y"^ < 1, the situation is very simple. The Cauchy 
potential u{z) outside fi, as was noted earlier, equals 



\z\ > 1, 



where c is a constant. Hence, outside the lens equation becomes 

c 

z — - 
z 



7z — w, 



(7.6) 



(7.7) 



a well-known Chang-Refsdal lens (cf., e.g., [1]) that may have at most 4 
solutions except for the degenerate case 7 = w = 0, when the Einstein ring 
appears. In particular, when 7 = 0, u; ^ 0, such mass distribution may only 
produce two bright images outside il. For z : [zj < 1 inside the lens the 
potential is still calculated by switching to polar coordinates: 



ip(r)rdrdO 







27T 



ip(r)dr 



rde 



re'- 



2tt 



ip(r)rdr 



d0 



re 



|z| 





27r 

z 



ip(r)rdr. 



In particular, for the "isothermal" density ip{r) 

u{z) = — \z\, 



so the lens equation (7.7 1 becomes 



z--^|z| 
z 



■ 7Z 



i, (pl yields for z : |z| < 1 



(7.9) 
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where c is a real constant. Equation (7.9 1 can have at most two solutions 
inside (only one, if 7 = 0), again, excluding the degenerate case of the 
Einstein ring. Furthermore, since Burke's theorem allows only an odd number 
of images, the total maximal number of images for an isothermal sphere 
cannot exceed 5 (4bright+l dim) as before (or < 3, i. e.,(2+l) if 7 = 0). Note, 
that strictly speaking, Burke's theorem cannot be applied to the isothermal 
density because of the singularity at the origin. Yet, since the density is radial 
and smooth everywhere excluding the origin and because it is clear from ( 7.9 ) 
that the origin cannot be a solution, Burke's theorem does apply yielding the 
above conclusion. 

The problem of estimating the maximal number of "dim" images inside the 
lens formed by a uniform mass-distribution inside a quadrature domain f2 
(cf. §4) of order n is challenging. In this case the Cauchy potential in (4.1) 
inside equals to the "analytic" part of the Schwarz function S{z). It is 
known that S{z) is an algebraic function of degree at most 2n. Yet, the sharp 
bounds, similar to those in Theorem |3.1| for the number of zeros of harmonic 
functions of the form z — a(z), where a{z) is an algebraic function, are not 
known. 

Another interesting and difficult problem would be to study the maximal 
number of images by a lens consisting of several elliptical mass distributions. 
Some rough estimates based on Bezout's theorem can be made by imitating 
the calculations in §5. Yet, even for 2 uniformly distributed masses these 
calculations give a rather large possible number of images (< 15) while, so 
far, only 5 images by a two galaxies lens and 6 images by a three galaxies 
lens have been observed - cf. [6,23]. 
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